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By relaxing the conventional assumption of a purely gravitational interaction between dark energy 
and dark matter, substantial alterations to the growth of cosmological structure can occur. In this 
work we focus on the homogeneous transfer of energy from a decaying form of dark energy. We 
present simple analytic solutions to the modified growth rates of matter fluctuations in these models, 
and demonstrate that neglecting physics within the dark sector may induce a significant bias in the 
inferred growth rate, potentially offering a false signature of modified gravity. 
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I. INTRODUCTION 

One of the major goals of future cosmological studies 
is to categorise dark energy as one of three candidates: 
a cosmological constant, a physical fluid, or simply some 
manifestation of new gravitational physics. Geometric 
measurements such as supernovac and baryon acoustic 
oscillations will scrutinise the null hypothesis of a cosmo- 
logical constant. However, a physical fluid and modified 
gravity could both reproduce almost arbitrary expansion 
histories, so distinguishing them requires a study of the 
growth of structure. It is the validity of this diagnostic 
step that provides the focus of the present work. 

Anisotropic stress in the dark energy fluid may mimic 
metric theories of gravity, as demonstrated by Kunz & 
Sapone [1] and more generally by Hu [2 . Here we take 
this one step further: even without anisotropic stress, 
it can be shown that significant deviations in structure 
growth are achievable. The only requirement is that some 
interaction should exist between dark matter and the 
smooth dark energy component, permitting the transfer 
of energy. 

Any change in the dark energy density is convention- 
ally attributed to the equation of state w{z), which dic- 
tates the adiabatic behaviour of a physical fluid. But den- 
sity evolution may also arise from non-gravitational inter- 
actions with other fluids, thereby violating adiabaticity. 
Cosmologies with energy exchange have been extensively 
studied in the literature [5l411j. and generate an expan- 
sion history that is fully reproducible by a single inert 
scalar field whose evolution matches the effective equa- 
tion of state Wf,s{z)- The degeneracy between these two 
models may be broken by studying the growth of struc- 
ture, which is disrupted by the evolving matter density 
[l2] . Yet it is this same test that would conventionally 
be used to identify modifications to gravity. 

In ^Hjwe begin by reviewing a simple case of dark en- 
ergy decaying into a form of dark matter. The evolution 
of perturbations are quantified in pil[ before exploring 
the observational consequences in pV[ 



II. ENERGY EXCHANGE 

When speculating on possible interactions amongst the 
lesser known constituents of our Universe, we are span- 
ning a remarkably broad class of models, with potentially 
numerous degrees of freedom. This could extend to new 
regimes of dark physics, such as dark matter particles 
spontaneously decaying into a relativistic dark species 
such as neutrinos or massless particles, e.g. the dark pho- 
tons speculated by Ackerman et al |13| . This particular 
example would be compatible with current observations 
provided the dark matter particle is sufficiently massive, 
and the dark fine-structure constant is sufficiently small. 
For the remainder of this work, we shall focus on the case 
of dark energy decaying into a form of dark matter, and 
explore the observational consequences. 

As an illustrative example, we study the simple case of 
a cosmological (almost) constant that decays into a ho- 
mogeneous form of dark matter. The selection w — —1 
for the dark energy fluid bypasses the various instabil- 
ity issues highlighted in previous work [5J [71 [5] . Instead 
we choose to focus on the behaviour of the dark matter 
perturbations. The evolution of the mean matter density 
Pm is dictated by the conservation equations 



Pa 



Pm + "y^Pm = Q, 



(1) 



(2) 



where the prime denotes a derivative with respect to con- 
formal time, and H = al ja = d is the conformal Hubble 
parameter. The interaction parameter Q controls the 
rate of energy transfer, and is generally considered to 
be a function of either the dark matter or dark energy 
density. 

The central result of this paper will be to demonstrate 
that models of interacting dark energy modify the growth 
of large scale structure, with an explicit example that 
exhibits a constant decrement c, such that 
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thus providing another mechanism for anomalous growth, 
aside from modified gravity or anisotropic stress. The 
magnitude of c is determined by both the nature and 
strength of the interaction. 

We adopt an empirical modification to the density evo- 
lution, which reproduces the functional form for a num- 
ber of models such as interacting quintessence [21 [H], 
given by 



Pm — PrnOO- 



-3+e 



(4) 



where the parameter e <C 1 dictates the rate of energy 
transfer. This particular parameterisation is advanta- 
geous in providing the necessary scaling of the dark en- 
ergy density to help resolve the naturalness and coin- 
cidence problems. The dark energy density maintains 
a magnitude comparable to the dark matter density, al- 
though it does not help explain why the transition to dark 
energy domination has occurred at recent times. Further- 
more, the simple form of Q allows significant progress 
to be made analytically. In this case the evolution of the 
dark energy density is, for a flat universe, given by 



A. Background Dynamics 

Here we generalise the derivation of linear structure 
growth by Linder & Cahn [15] to incorporate the differ- 
ent background evolution of both dark energy and dark 
matter, pm cx a"'^"'"'. Including energy exchange at the 
background level only is insufficient, as we will show, but 
it is instructive to identify the contribution this yields in 
the following derivation of /(a) = dlnJ/dlna. Starting 
from the differential equation 



dt2 



2i7(a) — -47rp,„J = 0, 



dt 



(8) 



which may be rewritten as 
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(9) 



(a)=0, 



and utilising 



Pa = PAQ + 



3-e 



PmO [o 



-3+e 



1] 



(5) 



At high redshifts, this leads to the tracking behaviour 
characterised by 



flAia -> 0) 



e 

3' 



(6) 



until we approach the era of radiation domination. Thus 
in effect this model adjusts the dark energy decay life- 
time with epoch; for a constant lifetime, we would expect 
Pa ^ constant at high z and pA at late times. Our 
analysis treats the parameter e as a constant, though re- 
mains valid provided e varies sufficiently slowly, satisfying 
the condition de/dlna <C e. 

Given the dark matter mass evolution Q , one can see a 
direct relation between our phenomenological approach, 
and that of a coupled scalar field as outlined in [5] . 
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III. GROWTH OF STRUCTURE 



Any form of non-gravitational interaction in the dark 
sector may be expected to impact upon the growth of 
structure. As we shall see, a decaying form of dark energy 
generates three distinct mechanisms for slowing the rate 
of structure growth. We quantify each of these in turn, 
starting with the modification to the expansion history 
Hiz). 



HVij2 ^ n^a-''+^ [1 + nA{a)/n^{a)] , (10) 
yields, to first order in (/(a) — 1) [15], 



where 



l--r!A(a) 



I{a) = -ia-(5+^)/2 r ^^(a')(5+0/2f]A(a') • (12) 

A I - n' 



(11) 
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Before evaluating this integral, we first need to consider 
the scaling behaviour of dark matter and dark energy 
at early times. To proceed, we simply recast the dark 
energy density ^ in the form 

nA{a) = ni+ n2(a) + O (e^) + O {nlia)) , (13) 

where fli is the constant component, and 5^2 is the col- 
lection of terms that follow a power law: 



0.2 = 



(14) 



provided we assume the contribution from radiation is 
negligible. 

Using these redshift dependencies in /(a), then to first 
order in deviations from matter domination we find: 



I{a) = - 



11 - 3e 



r!2(a) 



(15) 



3 



Writing J72(a) = ^A{a) — leaves us with 



I{a) 



30 66 



22 



+ 0{e')+0{nlia)) . (16) 



parameterisation of energy exchange, the density and ve- 
locity perturbation equations are greatly simplified from 
the general case presented in [5]: 



Substituting our solution back into (111 gives 



/(a)-l-(^-4e)llA(«)-i. (17) 

Of these two new terms involving e, the second is of 
greater importance. This contrasts with the conventional 
approximation for /(a) |161 117j given by 



/(«) 



1 - J^A{a) . 



(18) 



If we insist on maintaining the definition 7 « 
[1 — /(a)] /flAia), the extra ||e will mean 7 is no longer 
a constant to first order in llA(a). This suggests that a 
constant 7 is no longer a viable approximation for solving 
the perturbation equations. Rather than face the diffi- 
culty in generalising to a redshift dependent 7, we note 
that the constancy of 7 can be maintained at the ap- 
propriate level of approximation by simply changing the 
approximate solution such that 



g(a) =e^o"(d-7a')K(-')--i+b] 
The new constant b means that 

7«-(/(a)-l + fe)/r!A(a), 



(19) 



(20) 



and we see that if b is chosen to offset the constant terms 
in (/(a) — 1), then 7 will once again be a constant. For 
the case of pure background energy exchange, 



6' + eHS + 9- 3$' = 0; 



(24) 



e' + ne -k^"^ = 0. (25) 

The differential equation governing /(a) now becomes 



d/ /dlnff^ 
d In a 2 d In a 



/(/ + 2) + 2e 



2 d In a 2 
and carrying these new terms through gives 



(26) 
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1 - ^n^{a) 



(1 - 4e)/(a) . (27) 



The integral I (a) has already been evaluated (16), and 
so we find: 



or 



where 



23 

f{a) = ni{a)--e, 
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(30) 
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and hence our first modification to the growth rate is 



fia) c nZ^ia) - 



12e 

5^' 



where 7 is subject to a minor perturbation 
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C. Inertial Drag 

In the previous subsection we implicitly assumed that 
there is a pure energy transfer, from the perspective of 
the dark matter rest frame. If however the energy trans- 
fer occurs in the CMB frame, the newly formed dark mat- 
ter will not be instilled with the appropriate bulk motion, 
and a drag term appears in the perturbation equations. 

The density and velocity perturbation equations are 
then of the form 



6' + eHS + 9- 3$' = 0, 



(31) 



B. Dilution 

The energy exchange does not only affect the back- 
ground dynamics; it also generates new terms in the per- 
turbation equations, so ^ is no longer valid. By adding 
a homogeneous contribution to the background density, 
the fractional contrast, 6, is suppressed. In this particular 



6l' + -H(l + e)6'-fc2* = 0. 



(32) 



The physical interpretation of the eUS term in (31 ) is 



one of dilution, as addressed in the previous subsection. 
A second effect which also slows structure growth is the 
extra drag term in (32), eHO. This arises due to the 
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FIG. 1. The fractional change in the evolution of linear per- 
turbations in the presence of a decaying cosmological constant 
with e — 0.01 (upper) and e = 0.04 (lower), as defined in Q. 
The pair of dash-dotted lines highlight the slight suppression 
of growth induced by the diluting effect of transferring ho- 
mogeneous energy into the dark matter frame. The dashed 
lines correspond to an energy transfer in the CMB rest frame, 
with the extra deceleration arising due to the introduction of 
stationary matter. These are well described by the solid lines, 
which illustrate the analytic solution given by ( 38 1 . 
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FIG. 2. The logarithmic growth rate of both cold dark mat- 
ter (dashed), and the baryons (dotted), in the presence of 
a decaying cosmological constant with e — 0.04. The solid 



line illustrates the approximate solution given by (41 1, which 
traces the total matter perturbation. 



FIG. 3. The baryonic bias induced by e = 0.01 and e = 0.04. 



where we have neglected and the gravitational source 
term is rescaled as 



(35) 



These two new equations, ( 33 ) and ( 34 ), are precisely the 



equations for perturbations as if no energy exchange were 
taking place beyond a background level. The previous 
section provided us with an approximation solution in a 
universe with background energy exchange: 



After rescaling: 



or simply 



and we return to the solution (231 to find 
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7 = e . 
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(36) 



(37) 



(38) 



(39) 



appearance of stationary matter, reducing the mean flow 
rate. 

If we define new perturbations 5 = Sa'^ and = ^a*^, 
then substituting these rescaled variables into (31 1 and 
( 32 1 leaves us with 



~5' 



= 0, 



e' + 'He~kH ^Q. 



(33) 



(34) 



The physical mechanisms contributing to this new 
growth rate are summarised in Table [Tj The quoted total 
corresponds to the model where Q is unperturbed, and 
with the energy transfer occurring in the dark energy rest 
frame (that of the CMB). If instead the dark energy de- 
cay is sensitive to local fluctuations in the matter density, 
5Q oc T-LSprm the dilution term is no longer present. Sim- 
ilarly, if the energy exchange were to take place in the 
dark matter rest frame as opposed to the dark energy 
rest frame, the drag contribution vanishes. 
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Contribution 


A/A 


A7/e 


Background — 


/I 1 \ 
(5 + 55) 


D 

~55 


Dilution 


1 

~5 


10 

^55 


Drag 


4 

^5 


10 
+ 55 


Total 


67 

55 


6 

55 



TABLE I. Summary of the contributions to the growth 
rate from a decaying dark energy component, in units of e. 
Whether the effects of dilution and drag arise will depend on 
the behaviour and orientation of the interaction parameter 
Q. Note that the terms in the left column, those contributing 
directly to f{z), have a considerably greater impact than the 
corrections to 7 in the right column. 



D. Baryonic Correction 
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FIG. 4. The solid contours provide an example of the bias 
which may be induced in the gravitational growth index 7 
when making the false assumption that dark energy is stable. 
The true model, as indicated by the black dot, corresponds to 
e = 0.01. The dashed contours demonstrate the modification 
to the growth index induced by the elastic interaction model 
outlined in [TH] . 



Thus far we have assumed the universe comprises solely 
of dark matter and dark energy, but before considering 
the observational consequences we must first extend our 
analysis to include baryons. We therefore divide the mat- 
ter component into cold dark matter and baryonic terms, 
denoted pc and pb- The baryonic perturbations are sim- 
ply given by 



(40) 



The analytic solution in ( 38 ) now requires a minor 



rescaling to compensate for the introduction of baryonic 
matter. A simple weighting of the decay parameter e by 
the fraction of mass to which it applies is sufficient. 



67. 
55 



(41) 



(42) 



In Figure [2j this formalism is seen to provide an excel- 
lent description for the combined matter perturbation. 
Although we note that a more thorough treatment of the 
perturbations is required at very high redshifts {z > 100), 
where the radiation energy density starts to become sig- 
nificant. At later times the baryons closely track the 
growth rate of the dark matter perturbations, but with 
a slightly enhanced density contrast. Defining the bary- 
onic bias b = Sb/Sc, the present day value of this ratio is 
approximately given by 



as illustrated in Figure [3j 



(43) 



IV. OBSERVATIONAL CONSEQUENCES 

We now review their impact on cosmological observ- 
ables, using redshift-space distortions as a measure of 
the growth rate, before considering implications for the 
Integrated Sachs Wolfe (ISW) effect. 



A. Redshift Space Distortions 

In order to highlight the observational consequences of 
this energy exchange, we evaluate the appropriate Fisher 
matrix for a redshift survey ai z = 0.5 (see [19'), and 
combine this with the DETF Fisher matrix for Planck. 
We marginalise over the parameter set 



[wO, Wq, SIa, rifc, rife/l^, TT-s, As, /3, 7, dp 



(44) 



The standard cosmological parameters are taken to have 
fiducial values as derived from WMAP5 [5D]. In order 
to remain consistent with the CMB fisher matrix, when 
perturbing e we ensure the value of Pm at z = 1100 is held 
fixed. The Hubble parameter h is slightly perturbed for 
the purposes of distance estimation. We neglect the small 
change induced by the Integrated Sachs Wolfe effect (see 
TVB|. 



As a parameterisation of modified gravity, the growth 
of structure is taken to follow the form given by [1^1 [TT] 



(45) 



Clearly one ought to expect that, for the case of in- 
teracting models, applying the above prescription would 
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lead to a biased estimate of the growth index 7. In addi- 
tion to the extra e term in ( 38 1 , the value of 7 is biased 



if the evolution of ilm{z) does not run as expected. This 
shift is illustrated by the solid contours in Figure |4] which 
are generated with e — 0.04. 

In practice, we must deal with the Fingers of God in 
greater detail than a single parameter cTp. It should also 
be noted that these models may have some impact on 
virialised structures, though we leave this as a topic for 
future investigation. 



B. Integrated Sachs Wolfe Effect 

One restriction that dark energy models must satisfy 
is not to overpredict the Integrated Sachs Wolfe effect. 
Excessive change in the gravitational potential would in- 
variably generate large fluctuations in both the CMB and 
CMB-LSS cross-correlation on large angular scales. In- 
deed this has already been employed by Valiviita et al. [B] 
to assist in constraining models of decaying dark energy. 

Since we are still working in the context of General 
Relativity, the gravitational potential is readily evaluated 
via the Poisson equation 



$(fc,a) = -47rGp(t)a' 



2-5(0) 
A:2 



(46) 



(47) 



This evolution is illustrated in Figure [5] for the standard 
case of flat ACDM, alongside small perturbations in the 
decay parameter (e) and global curvature (fJfc). While 
there is a distinctive change in behaviour at high red- 
shift, this is a regime that lies out of reach for cross- 
correlation studies. Radiation may also start to become 
significant at this point, we have not included this effect 
in our treatment. 

The ISW signal is, at best, a ~ 5cr observation, and as 
such there is significant room for flexibility in the antici- 
pated signal strength. Taking a 20% change as the max- 
imum permissible, this corresponds to an approximate 
upper limit of e < 0.1. This appears broadly consistent 
with the findings of Valiviita et al. [6^ , who established an 
upper bound of |r| < 0.23i/o for a constant dark energy 
decay rate, from a combination of WMAP, supernovae, 
and BAO data. The two parameterisations are related 
by r = —"He. Our earlier analysis on the growth rate 
may be naturally extended to this model, which we also 
find to be well described by the prescription 



A/ cx T/H . 



(48) 
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FIG. 5. The decay rate of the gravitational potential as expe- 
rienced in a flat ACDM Universe (sohd), and with perturbed 
cosmologies e = 0.01 (dashed) and 57^ = —0.01 (dotted). 



V. DISCUSSION 

In the case of a decaying cosmological 'constant', we 
have demonstrated that the growth rate of large scale 
structure is subject to a constant decrement / f77„ — c. 
For the simple model we consider, this can largely be at- 
tributed to the dragging effect on bulk motions induced 
by the production of stationary matter. Smaller contri- 
butions arise from both the background dynamics, and 
the diluting effect of gradually introducing a homoge- 
neous density field. It is also important to note that 
in more general interacting models the presence of dark 
energy perturbations may also influence the growth of 
structure. These act to enhance the growth rate, and 
may overpower the mechanisms considered here. Indeed 
in many cases the dark energy perturbations in coupled 
models lead to pathological instabilities. 

Energy exchange within the dark sector leads to a 
change in the comoving matter density, which in turn in- 
vokes a number of changes to cosmological observations. 
If these changes are not taken into consideration, then 
suppression to the growth of large scale structure leads to 
a naively inferred value of the growth index 7 rising above 
the conventional value 7 > 0.55. While the physical mo- 
tivation for such models remains unclear, we believe this 
is no less true for current approaches to modified gravity. 

As noted by Blandford et al. [21], a blind cosmolo- 
gist living in the radiation dominated era might measure 
the evolution of the scale factor and erroneously con- 
clude that the expansion is driven by a scalar field with 
an exponential potential. This thought experiment may 
be extended further: a "dark" astronomer only capable 
of studying the dark matter would notice a strange be- 
haviour of non-linear structure. They might attempt to 
construct a modified theory of gravity which accounts for 
this behaviour. However the true source of this discrep- 
ancy is the momentum exchanged between the baryons 
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and photons, inducing new features to the growth of 
structure. This modification is apparent both at the era 
of recombination, and the present day. We are at risk of 
falhng into a similar trap here: the baryocentric assump- 
tion that neither dark energy nor dark matter exhibit 



any complex behaviour may lead to dark physics being 

mistaken for modified gravity. 
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